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■ Abstract 

| Integral expressions are determined for the elastic displacement and stress fields due to 

I ■ stationary or moving dislocation loops in finite samples. These general expressions are valid 

for anisotropic media as well. Specifically for the stress fields, a line integral representation is 
t-H , found, thus showing rigorously the independence of the stress fields with respect to the choice 

^ ■ of slip planes. In the stationary case the line integral representation involves calculating 

CO ' a "vector potential" dependent on the specific geometry of the sample. Two examples of 

geometries, isotropic half space and thin plate, are shown where the "vector potential" has 
been explicitly determined. With this general method one recovers some earlier specific results 
in these geometries. 

t> : 

1 Introduction 
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S A general formula for the displacement and stress fields generated by a dislocation loop undergoing 
arbitrary motion in an infinite medium was obtained some years ago by Mura (1963). It is written 
as a convolution of the medium's impulse response with a source localized along the dislocation 
loop (i.e. independent of the loop's slip plane). These formulae do not apply to the case of finite 
samples since they are derived assuming homogeneity in space, i.e. a dependence of the Green's 
function in x — of, which indeed it is no longer true in a finite sample. The purpose of this paper 
is to announce a generalization of these formulae to finite samples, as well as explicit forms for the 
geometries of a half space and thin plate. The main results will be presented, with some examples, 
5_i ■ but the detailed proofs will be presented elsewhere. 

Motivation for this work comes partly from an ongoing project that attempts to understand 
recently observed dynamic instabilities for cracks in thin plates (Sharon, Gross and Fineberg 1996, 
Boudet, Ciliberto and Steinberg 1996) in terms of the interaction of the plate's oscillations with 
the crack's tip regarded as a continuous distribution of infinitesimally small dislocations (Lund 
1996). An intuitive thought behind this work was that the stress fields produced by stationary or 
moving dislocation loops in a finite sample should not depend on the slip planes chosen, and that 
indeed the stresses will be continuous through the slip planes, as it was in an infinite medium. 
This is rigorously shown here: the expressions found for the stress fields involve only line integrals 
over the dislocation loops. 

When the loops are stationary in order to find the explicit form of the line integrals appearing 
in the stress fields, one should calculate a "vector potential" appropriate to the finite sample in 
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question. This was done for the two examples shown here: a half space and a thin plate. In a 
non stationary situation the displacement velocity field can be written in terms of line integrals, 
and from this one deduces that the temporal partial derivative of the stress field can be written 
in terms of them too. One finally deduces that the time dependent stress field can be written in 
terms of line integrals since one knows that the initial condition (stationary case) and its temporal 
derivative can both be written in terms of line integrals. 



2 Displacement fields due to dislocations in finite samples 
2.1 Stationary case 

The following formula is valid in a finite sample for the displacement field produced by a stationary 
dislocation loop (derived in detail elsewhere). It is written as an integration over the slip plane 5", 
and involves a discontinuity of the displacement there, given by the Burgers vector bf. 

U m (x) = -bi f dS'ja^ix^x). (1) 

JS' 

a^ m is the elastic stress associated with the static Green's function Gi m (x, x'), which is the dis- 
placement (direction (I)) produced by a localized impulse force in space (x r ) (direction (to)): 

cm dS^ Gimix,x ' ] = i^ m{x ' x,) = ~ SimS{x x>) ■ (2) 

It satisfies the free surface boundary condition: 

<jC m (xs,x')n 1 (x s ) = 0, (3) 
with rij(xs) the components of the normal to the free surface at its point xs- 



2.2 Non stationary case 

The following formula is valid in a finite sample for the displacement field produced by a moving 
dislocation loop. It involves integration over past slip planes S(t'): 

U m (x,t) = -bi[ dt' [ dS'^f^ix' ,x;t-t') . (4) 

J-oo JS(t') 

The Green's function Gi m (x, x';t — t') is the displacement (in direction (i)) produced by a localized 
impulse force (in direction (m)) in space (x') and time (t 1 ): 

P Q -^T^ x 'i * " " ^T* 7 ^' x '> *-*') = S im 6{x - x')8(t - t') , (5) 
and satisfies that the normal stresses it produces at the free surface of the sample are zero. 
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3 Stress fields due to dislocation loops in finite samples 
3.1 Stationary case 

Starting from Eq. ([!]), the stress field produced by the static dislocation loop is: 

a pq (x) = -hC vqlm J dS'—o^tf, x) . (6) 

This can be expressed as a line integral along the dislocation loop: 

a pq (x) = hC pqlm [ dl' s {e sjl a% m (x',x) - A^ m (x',x)} , (7) 
v 



with the "vector potential" A\ (x',x) to be determined for each specific sample. It is defined by: 

^"V'Z) + ^°«V.*) = ^rs^ m {x',x) . (8) 

3.2 Non stationary case 

Since the problem is homogeneous under time translations it is possible to obtain an expression for 
particle velocity, that involves an integral along the dislocation loop only. Indeed, from Eq. ([!]): 

^(x,t) = b i f°° dt'Af dS'^^x'^t-t')}- f dS^ m (x\x;t-t')}. (9) 

The first term is zero since the Green's function vanishes at t — t' = ±oo. Since f dS ( t iy dt , dS'j = 
£ jpq IlW) dl' q V p (x' ,t') {L{t') is the dislocation loop bounding the slip plane S(t'), V p (x',t') is the 
loop's local velocity and ej pq the completely antisymmetric tensor in three dimensions), we have: 

CjJJ pOO f 

' m -{x.t) = -b l dt 1 dl'<T% m (x',x;t-t')e jpq V p (x',t'). (10) 



dt 



L(t') 



From this expression for the displacement velocity field, one obtains an expression for the partial 
time derivative of the stress field that involves line integrals over the dislocation loop: 

JU P9 (:M) = -C pqkm bi J dt' J ^ dl' q -^-a^ m {x',x;t-t')e jpq V p (x',t') . (11) 

As mentioned in the Introduction, this result together with the result for a static loop (consid- 
ered as an initial condition) mean that the time dependent stress field can be written in terms of 
integration of line integrals over the dislocation loop, consequently it is independent of the choice 
of slip planes. 

4 Examples for the stationary case 
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4.1 A screw dislocation in a half space 

Using the method in Eshelby and Stroh (1951) the displacement and stress fields in polar coordi- 
nates due to a screw dislocation lying along the z axis perpendicular to the free surface of a half 
space, can be obtained as (the half space exists for z > 0): 

u, m = y 

a z ,{r,z) = !°° dke-<".h{kr) 

Zirr Z7T J 

<J r e(r,z) = f dker kz J 2 {kr) (12) 



2tt ... 

(Ji{z) the Bessel's function of order one). Applying the general formula of Eq 



a zv (R, z) = bfie ari I dz'<rf*{z', z; -R) - bfi j dR' Af s (z' = 0, z; R' — R) , (13) 



with greek indices indicating coordinates x or y only, R = xi + yj, fi the shear elastic constant 

A Zt 



and C is an arbitrary curve on the free surface. The "vector potential" Ag zS at the free surface is 



given in Fourier space by: 

Af 5 {z' = 0, z; k) = ^-kpe Sri k v . (14) 

If one calculates a z $ (r, z) using Eq. ( |l3| ) , one reproduces the first and second terms in Eq. ( [l2| ) 
with the equivalent terms in Eq. (|l3|). Gosling and Willis (1994) have given a different line integral 
representation specifically determined for the geometry of an isotropic half space. 

4.2 A screw dislocation in a plate 

In Eshelby et al. (1951) the exact displacement and stress fields in polar coordinates due to a 
screw dislocation along the z axis perpendicular to the free surfaces of a plate, were obtained as 
(the plate exists for — h < z < h): 

1° 



U z (9) 



b f dk sinh(fcz) 
Ue{r,z) = - — / Ji(fcr) 



27T J k cosh(kh) 
fib fib f°° cosh(fcz) 
2-kt 2n J cosh(kh) 

^ e(r ' z) = 2^X *c^h(fo) J2(M • (15) 
The calculation of cr z o(r, z) using Eq. (0) is analogous to the half-space case. One needs the "vector 

zz ~ 





potential" A% zS evaluated at the surface, which in Fourier space is 



A»*S(+h -y-k\- k ^s,K (+ cosh(kz) sinhjfez) 
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